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Abstract 

Context. The formation and evolution of disk galaxies are long standing questions in Astronomy. Understanding the properties of 
globular cluster systems can lead to important insights on the evolution of its host galaxy. 

Aims. We aim to obtain the stellar population parameters - age and metallicity - of a sample of M31 and Galactic globular clusters. 
Studying their globular cluster systems is an important step towards understanding their formation and evolution in a complete way. 
Methods. Our analysis employs a modern pixel-to-pixel spectral fitting technique to fit observed integrated spectra to updated stellar 
population models. By comparing observations to models we obtain the ages and metallicities of their stellar populations. We apply 
this technique to a sample of 38 globular clusters in M31 and to 41 Galactic globular clusters, used as a control sample. 
Results. Our sample of M31 globular clusters spans ages from 150 Myr to the age of the Universe. Metallicities [Fe/H] range from 
-2.2 dex to the solar value. The age-metallicity relation obtained can be described as having two components: an old population with 
a flat age-[Fe/H] relation, possibly associated with the halo and/or bulge, and a second one with a roughly linear relation between age 
and metallicity, higher metallicities corresponding to younger ages, possibly associated with the M31 disk. While we recover the very 
well known Galactic GC metallicity bimodality, our own analysis of M3 1 's metallicity distribution function (MDF) suggests that both 
GC systems cover basically the same [Fe/H] range yet M31's MDF is not clearly bimodal. These results suggest that both galaxies 
experienced different star formation and accretion histories. 

Key words. Galaxy: globular clusters: general - Galaxies: star clusters: M31 



1. Introduction 

Globular clusters (GCs) are widely considered as excellent as- 
trophysical laboratories. Their age and metallicities, in particu- 
lar, trace the main (astro )physical processes responsible for the 
formation and evolution of their host galaxies. However, differ- 
ently from the Milky Way, most of the extragalactic GCs cannot 
be resolved into individual stars. 

The techniques to estimate ages and metallicities of ex- 
tragalactic GCs have classically fallen int o two broad cate gor- 
ies. Those relying on photometry (e.g. iFan et al 1 12006h are 
more susceptible to the age-metallicity degeneracy in the sense 
that young metal-rich populations are photometrically indis- 
tinguishable fr om older metal-poor populations (but see, e.g., 
iMa et al.ll2007l) . The different spectroscopic methods, however, 
are inspir ed by the Lick/IDS system of absorption line in- 
dices (e.g. lWorfhev et ai1ll994l and references therein) — fr om 
linear metallicity calibrations (e.g. iBrodie & Huchralll990h to 
those which perform the simultaneous ^-minimisation of a 
large number of spectral indices (see, e.g.. iProctor et al.ll2004h . 
Nevertheless, both photometric and spectroscopic methods are 
dependent on the acc urate modelling of simple stellar popula- 
Bruzual & Chariot 2003; Le Borgne etalj 
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Maraston 2005 
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ICid Fernandes & Gonzalez Delgadol (1201 ol) and references 
therein. This method has advantages such as making use of all 
the information available in a spectrum (making it possible to 
perform analysis at lower S/N) and not being limited by the 
physical broadening, since the internal kinematics is determined 
simultaneously with the population parameters. In some of its 
flavours, this method is also insensitive to extinction or flux 
calibrating errors. It has also been shown that full spectrum 
fitting reproduces better the results fro m colour-magnit ude 
diagrams (CMDs) than other methods (e.g. lWolf et al.ll2007l) . 

In this study, we have used the spectrum fitting code t/Ly&sQ 
dKoleva et al.l2009l) to compare, on a pixel-by-pixel basis, the in- 
tegrated spectrum of 38 spectra of GCs in M31 to simple stellar 
population models in order to derive their ages and metallicities. 
O ur sample comprises 3 5 integrated spectra previously analysed 
by Beasle v et all d2004l). and three outer halo clusters taken from 
lAlves-Brito et al.l d2009l) . In addition, we have al so analysed in- 
tegrate d spectra of 41 Galactic GCs presented in lSchiavon et al.l 
(120051) . whose populations have been studied with CMDs and 
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spectroscopy of individual stars. Both, Galactic and M31 GCs 
were analysed in the same way. The Galactic sample acts as a 
control sample to estimate the reliability of the fitting method. 

Our nearest (~780 kpc, iHollandl 1 1 9981) giant spiral galaxy 
M3 1 is a natural target to test the current ideas about the forma- 
tion a nd evolution of ga l axies in the Local Universe. As pointed 
out in lAlves-Brito et al I d2009l a remarkable difference between 
the M31 and the Galactic GC system is that M31 hosts more 
GCs (by a factor of 2-3) than the Milky Way. Furthermore, there 
have been suggestions that M3 1 may have a significant popula- 
tion of young and intermediate age (less than 8 Gyr) that is not 
found in the Galactic system. The metallicity distribution func- 
tion (MDF) of M 31 (through both giant stars and GCs) and o ther 
galaxies (see, e.g. lAlves-Brito et"aUl201 lfclUsher et al.l2012l) has 
been a controversial topic regarding its shape and distribution 
when compared with our own Galaxy. In particular, the pres- 
ence (or not) of colour- [Fe/H] bimodality is one of the interest- 
ing and intriguing open questions in the field. Therefore, these 
(dis)similarities between the different spiral galaxies in the Local 
Group need to be investigated through different techniques to 
better understand the different process(es) in which spiral galax- 
ies are formed and have evolved in the Universe. 

The paper is organised as follows. In Sect. 2 we describe our 
sample. In Sect. 3 we present the stellar population analysis. In 
Sect. 4 we discuss the results obtained. Concluding remarks are 
given in Sect. 5. 



2. Sample 

2.1. M31 GCs 

The in tegrated spectra used in t his s tudy were previously stud- 
ied by lAlves-Brito et alJ ri2009h and iBeaslev etaP d2004l) . The 
former provided spectra for three outer halo GCs at projected 
distances beyond « 80 kpc from M31. These spectra were ob- 
served with the cross-dispersed, high-resolution spectrograph 
HIRES instrument on the Keck 1 telescope, covering a broad 
wavelength range of AA = 4020-8520 A, at a spectral resolving 
power of R w 50, 000. Spectroscopic age and metallicities were 
obtained by using metallicity calibrations from Mgb, CH and 
Mg2 indices. In addition, the authors also employed the simul- 
taneous ^-minimisatio n of a large number of spectral indices, 
approach introduced bv lProctor et al.l d2004l) . 

The second source of data is detailed in Be aslev et al.l(l2004l) . 
who made a analysis of high-quality integrated spectral indices 
in M3 1 . We have only used the spectra obtained with the Low 
Resolution Imaging Spectrograph mounted on the Keck I tele- 
scope, covering a spectral range of AA = 3670-6200 A and with 
a full width at half-maximum (FWHM) resolution of 5 A. The 
sample chemical properties were studied through the measure- 
ment of Lick line strengths. 

In Table Q] we list the M3 1 GCs analysed in t he pr esent 
work, and we refer the reader to lAlves-Brito et al.l d2009l) and 
IBeaslev et al.l d2004l) for additional information about the ob- 
servations and data reduction. The projected distances shown 
in Table Q] were calculated relative to an adopted M3 1 central 
position of aj2ooo = 00 h 42 m 44 s .30, <5 J2 ooo = +41°16' 09.90"4. In 
addition, a position angle for the X-coordinate of 38° ([Kent et al.l 
119891) as well as a distance of 780 kpc dHollandlll99'8¥ were ad- 
opted. At this distance, 1 arc-min corresponds to 228 pc. 



2.2. Galactic GCs 

The sample of Galactic GCs spectra is given in Tab le |2l whose 
observations were taken from lSchiavon et al] d2005l) . The obser- 
vations were performed with a long-slit spectrograph in drift- 
scan mode in order to integrate the population within one core 
radius. The spectra cover the range AA = 3350 - 6430 A at a 
resolution of about FWHM = 3 A. The mean S/N varies from 50 
to 240 A 1 depending on the wavelength. We refer the reader to 
ISchiavon et al.l d2005h for more details on the observations and 
data reduction. 

The aim of analysing this Galactic sample is twofold: first, 
we use it as a control sample where we can evaluate the perform- 
ance of the fitting technique by comparing our results to stud- 
ies of CMD and spectroscopy of individual stars; and second, 
by ensuring that the same analysis is applied consistently to the 
M3 1 and Galactic sample, the latter acts as a reference system to 
which the M3 1 GC system is compared. 

We searched the literature for good independent determina- 
tions of ages and metallicities of this sample. Fitting theoretical 
isochrones to GC CMDs is generally accepted as the most se- 
cure age determination possible when using photometry; how- 
ever, the results do vary between sets of isochrones and meth- 
ods of analysis, and the absolute derived age depends on model 
zero points, input physics, colour- T e jf transfor mation, distance 
uncer t ainties and foreground reddening (e.g. IChabover et aP 
119981; iBuonanno et"aL1 Il998t iMeissner & Weissl 120061) In this 
work we relied on the relative a ges derived homogeneo usly 
for a large sample of clus ters from lDe Angeli etaP d2005l) and 
M ann-Franch et al. Those relative ages were converted 

to absolute ages by adopting a conserv ative value of 13 + 2.5 
Gyr for 47 Tucanae dZoccali et al.ll200ll) . We additionally added 
val ues forNGC 6528 and NGC 6553 from lMomanfetail d200l 
and lBeaulieu et al.l d2001l) . respectively. 

Regarding the metallicities, we adopt two homogen- 
eous compilations found in literature: the compilation by 
ICarretta et al.l (12009b . who analysed high resolution stellar spec- 
tra in 19 globular clusters, and brought the lower resolutio n 
measurements f rom Zinn & West! dl984l) : lKraft & Ivansl d2003l) : 
iRutledge et aP dl997^ to a common metallicity-scale; and the 
com pilation bv [S chiavon et al.1 ( 2005[). based on mea surements 
from lKraft & Ivansl d2003l) : ICarretta & Grattonl dl997l) . 



3. Analysis 

We obtained ages and metallicities for the GCs through the com- 
parison of their i ntegrated spectra to SSP models, using the pub- 
lic code ULySS (Koleva et al. 2009), described briefly below. 

ULySS is a software package performing spectral fitting in 
two astrophysical contexts: the determination of stellar atmo- 
spheric parameters and the study of the star formation and chem- 
ical enrichment history of galaxies. In ULySS, an observed spec- 
trum is fitted against a model (expressed as a linear combination 
of components) through a non-linear least-squares minimisation. 
In the case of our study, the components are SSP m odels. 

We adopt SSP models bv lVazdekis et al.l (l2010t) . which cover 
the wavelength range 3540 -7400 A at a resolution of FWHM 
~ 2. 5 A. Ages range between 63Myr and 18 Gyr, and metal- 
licities [Fe/H] between -2.32 and + 0.22 dex. These models 
are based on MILES st ellar library dSanchez-Blazauez et al.l 
120061: ICenarro et al.ll2007l) and lGirardi et al.l (l2000h evolutionary 
tracks. Being based on an empirical library, the models follow 
the abundance pattern of the solar neighbourhood, namely they 
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Table 1. M31 GCs analysed in the present work. 



ID 


R.A. (J2000.0) 


Decl. (J2000.0) 


d P 
(kpc) 


S/N 
(pixel -1 ) 


Reference 


B126 


00:42:43.481 


+41:12:42.47 


0.79 


73 


a 


B134 


00:42:51.678 


+41:44:03.42 


0.52 


70 


a 


B158 


00:43:14.406 


+41:47:21.28 


2.40 


120 


a 


B163 


00:43:17.640 


+41:27:44.91 


3.01 


170 


a 


B222 


00:44:25.380 


+41:14:11.62 


4.37 


30 


a 


B225 


00:44:29.560 


+41:21:35.27 


4.69 


169 


a 


B234 


00:44:46.375 


+41:29:17.77 


6.03 


71 


a 


B292 


00:36:16.666 


+40:58:26.58 


17.15 


32 


a 


B301 


00:38:21.581 


+40:03:37.16 


20.09 


31 


a 


B302 


00:38:33.500 


+41:20:52.29 


10.81 


28 


a 


B304 


00:38:56.940 


+41:10:28.41 


9.85 


34 


a 


B305 


00:38:58.922 


+40:16:31.78 


16.74 


15 


a 


B307 


00:39:18.477 


+40:32:58.05 


13.27 


20 


a 


B310 


00:39:25.752 


+41:23:33.14 


8.67 


35 


a 


B313 


00:39:44.599 


+40:52:55.05 


9.38 


52 


a 


B314 


00:39:44.599 


+40:14:07.95 


16.15 


17 


a 


B316 


00:39:53.604 


+40:41:39.29 


10.78 


20 


a 


B321 


00:40:15.545 


+40:27:46.50 


12.78 


30 


a 


B322 


00:40:17.270 


+40:39:04.70 


10.57 


69 


a 


B324 


00:40:20.477 


+41:40:49.38 


8.33 


71 


a 


B327 


00:40:24.107 


+40:36:22.52 


10.91 


30 


a 


B328 


00:40:24.529 


+41:40:23.15 


8.13 


30 


a 


B331 


00:40:26.642 


+41:42:04.28 


8.35 


15 


a 


B337 


00:40:48.477 


+42:12:11.06 


13.71 


178 


a 


B347 


00:42:22.892 


+41:24:27.58 


8.79 


61 


a 


B350 


00:42:28.442 


+40:24:51.12 


11.73 


39 


a 


B354 


00:42:47.645 


+42:00:24.73 


10.10 


30 


a 


B365 


00:44:36.445 


+42:17:20.57 


14.77 


60 


a 


B380 


00:46:06.238 


+42:00:53.13 


13.36 


43 


a 


B383 


00:46:11.498 


+41:19:41.48 


8.94 


85 


a 


B393 


00:47:01.204 


+41:24:66.39 


11.16 


43 


a 


B398 


00:47:57.786 


+41:48:45.66 


15.32 


44 


a 


B401 


00:48:08.508 


+41:40:41.93 


14.96 


54 


a 


MGC1 


00:50:42.459 


+32:54:58.78 


116.05 


20 


b 


MCGC5 


00:35:59.700 


+35:11:03.60 


78.49 


20 


b 


MCGC10 


01:07:26.318 


+35:46:48.41 


99.85 


20 


b 


NB16 


00:42:33.094 


+41:20:16.41 


1.05 


142 


a 


NB89 


00:42:44.780 


+41:14:44.20 


0.34 


110 


a 



Notes. — The first column gives the GC identification, as in the Bologna catalogue (Galleti et ai1 l2004h . Second 
and third columns give the coordinates of the objects, in J2000. Fourth column shows the projected distance of 
the cluster to the centre of M31 (see text in §2.1). Fifth column shows the S/N of the spectra as given in the 
corresponding paper; when this information was not available (clusters B302, B305, B307, B316, B331, B354), we 
estima ted the S/N around 5 000 A usin g standard IRAlQ routine s. Sixth column gives the reference of the spectra: 
(a) for lBeaslev etail d2004l) and (b) for lAlves-Brito et all d2009l) . 

" IRAF is distributed by the National Optical Astronomy Observatory, which is operated by the Association of Universities for 
Research in Astronomy, Inc., under cooperative agreement with the National Science Foundation. 



are solar-scaled around solar metall icities, and ff-enhanced for 
low metallicities dMilone et al.ll201ll) . 

ULySS matches model and observation continua through a 
multiplicative polynomial, determined during the fitting process. 
Therefore, ULySS is not sensitive to flux calibration, galactic ex- 
tinction, or any other cause affecting the shape of the spectrum. 
We run ULySS with its global minimisation option, i.e., each fit- 
ting was performed starting from several guesses; this is an im- 
portant feature that minimises the risks of results being biased 
by local minima. 



For each spectrum we run 300 Monte Carlo (MC) simula- 
tions. The simulations consist of analyses of the spectrum with 
a random noise added, accord i ng to the S/N of the observation. 
As explained in iKoleva et al.l (l2009h . the MC simulations take 
into account the correlation between pixels and thus reproduce 
the correct noise spectrum, giving a robust estimate of the er- 
rors. We adopt the final stellar populations parameters to be the 
mean values of the 300 simulations, and our final uncertainties 
are the 1 <x value of the simulations results. We quote a minimum 
error of 10% in age to account for possible systematic errors in 
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the isochrones used in the SSP models (e.g. lGallart et al.|[2005l) . 
when the 1 <x value is smaller than this limit. 

A possible disadvantage of spectral fitting, compared to Lick 
indices, is that the method may in principle be sensitive to the 
wavelength ra nge fitted. Reports o n this effect in literature are 
contradictory: iKoleva et al.l d2008l) has concluded that the sens- 
itivity to th e wavelength range i s not critical; at odds with this 
conclusion. IWalcher et al.l d2009l) reports that while fitting glob- 
ular clusters of the Galactic bulge, the results are more accurate 
when the fitting is limited to a wavelength range ~ 530 A wide 
centred at 5 100 A. 

We thus performed some tests on the dependence with the 
wavelength window by fitting the Galactic GCs sample at differ- 
ent wavelength ranges. The widest range we tested was 3650 - 
6150 A (the total range covered bv lBeaslev et alj|2004l observa- 
tions) and the shortest 4828 - 5364 A (favou red bvlWalcher et al] 
I2009h . We also tested the wavelength used in lKoleva et al.l (I2008L 
4000 - 5700 A) and the range 400 - 5400 A, similar to the cov- 
erage of lAlves-Brito et all (120091) observations. 

4. Results 

4.1. Stellar populations in the Galaxy — verifying the method 

By analysing the Galactic sample at different wavelength win- 
dows, we verified that the wavelength choice has little influ- 
ence on the metallicities derived via spectral fitting. On the other 
hand, the ages may change in a non-negligible way: in particu- 
lar, fitting the widest range resulted in nearly a third of the GCs 
being fitted with intermediate ages (down to 4Gyr). Possible 
sources of errors are: blue horizontal branches or blue stragglers 
not properly taken into account on the models, which is known to 
affe ct spectroscopic ages both in spectra l fitting and Lick ind ices 
(e.g. lde Freitas Pacheco & Barbuvll99llKoleva et al.l2008l) . de- 
ficiencies in the observations (diffuse light affecting the blue re- 
gion and po or subtraction of te lluric lines affecting the red, as 
mentioned in IKoleva et al.l2008l) and limitations in the SSP mod- 
els (chemical patterns different than the ones in MILES library). 

In the present work we favour the wavelength range 4000 - 
5400A which, among the wavelength windows we tested, is the 
one which best reproduced the CMD ages of the Galactic glob- 
ular clusters. The results of this fitting run are shown in Table 
|2]together with literature values. Comparisons between our res- 
ults and literature are also shown in Figures Q] and [2] for ages 
and metallicities, respectively. Several of the Galactic GCs have 
multiple observations and we fitted each spectrum individually. 

In Fig. Q] it can be seen that six clusters have spectroscopic 
ages which are outside the range of allowed CMD values: NGC 
2808 (GC #5 in Fig. 1), NGC 5286 (GC #7), NGC 6121 (GC 
#12), NGC 6388 (GC #23), NGC 6441 (GC #24), NGC 7078 
(GC #33). 

In the case of NGC 5286 (GC #7), it is striking that three 
observations of the same cluster resulted in ages different by al- 
most 6 Gyr. According to the observation log in lSchiavon et alj 
(12001 . the three observations of this GC were taken at differ- 
ent slit positions or extraction aperture. A visual inspection of 
the fitting showed that in two exposures the observed Balmer 
lines are visibly narrower than the model, even though the metal- 
lic lines are well fitted (this pattern is also present in NGC 
6752, GC # 32). We suspect that contamination from foreground 
stars hampered the observations. A visual inspection of the field 
around this GC in Aladii^ shows a handful of bright foreground 

2 http://aladin.u-strasbg.fr 



20 




I .... i .... i .... i .... i .... i .... i ... ■ 

5 10 15 20 25 30 35 
Cluster § 



Figure 1. Comparison between the ages derived via spectral fit- 
ting (filled circles) and ages derived via CMD (filled grey pat- 
tern) for our sample of Galactic GCs. The number in the x-axis 
identifies the GCs (1: NGC 104; 2: NGC 1851; 3: NGC 1904; 4: NGC 2298; 
5: NGC 2808; 6: NGC 3201; 7: NGC 5286; 8: NGC 5904; 9: NGC 5927; 10: 
NGC 5946; 11: NGC 5986; 12: NGC 6121; 13: NGC 6171; 14: NGC 6218; 15: 
NGC 6235; 16: NGC 6254; 17: NGC 6266; 18: NGC 6284; 19: NGC 6304; 20: 
NGC 6342; 21: NGC 6352; 22: NGC 6362; 23: NGC 6388; 24: NGC 6441; 25: 
NGC 6528; 26: NGC 6544; 27: NGC 6553; 28: NGC 6624; 29: NGC 6637; 30: 
NGC 6652; 31: NGC 6723; 32: NGC 6752; 33: NGC 7078; 34: NGC 7089). 



stars close to the cluster. Given that the observations were per- 
formed in drift-scan mode, diffuse light from one of these stars 
might explain the failure of fitting the line profiles. For the re- 
maining clusters, there was no obvious pattern in the visual in- 
spection of the fitting, so we further investigated other sources 
of problems that could impact the derived ages. 

It is well established in literature (e.g. iGratton et al.l 

12011 that virtually all clusters harbour populations with anti- 
correlated C-N, O-Na (and sometimes Mg-Al) abundances. In 
this sense, a comparison to strict SSP models migh t not be ad- 
equate in many cases. In ICoelho et al.l J20 1 lll20 1 21) the authors 
investigate, by means of stellar population modelling, the ef- 
fect of populations with CNONa variations in spectral indices 
and integrated spectra. Inside the w avelength range fitte d in 
this work, we obtain from Fig. 5 in ICoelho et al.l (1201 2l) that 
the range -4150-4220A is affected by CNONa anti-correlated 
abundances. We investigated if these chemical variations could 
explain the age mismatch by repeating the fitting of the clusters 
NGC 2808, NGC 6121, NGC 6388, NGC 6441, NGC 7078, 
masking the CNONa affected region. We verified that the mask- 
ing does not improve the age results, the differences between 
masking or not the CNONa region being Aage = 0.1 ± 0.2 Gyr. 
We conclude therefore that CNONa anti-correlations in the in- 
tegrated spectra of clusters cannot explain the cases where spec- 
troscopic ages do not match the CMD values. 

The test above does not guarantee against the effect of mul- 
tiple main sequences and/or sub-giant branches, clusters for 
which remarkably a single age isochro ne is not adequate to fit the 
CMD. This is the case of NGC 2808 dPiotto et al.ll2007l) . which 
is known to harbour three main sequences and a very complex 
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horizontal branch. In such striking cases, it is not disquieting that 
comparing the observations to SSP models would fail. 

Except for NGC 2808, the other clusters with deviating 
spectroscopic ages are all younger than the CMD ages by 
at least 2-3 Gyr. Our initial interpretation was the poten- 
tial presence of HB morphologies which are not well rep- 
resented in the SSP models. Extended HB morphologies are 



lower values (e.g. de Freitas Pacheco & Be 


rbuvl 


1995b Lee et alJ 


20001 Schiavonetal J 120041 Mendel et al.l 


2007 


: lOcvirkl 20101). 


Koleva et al. ( 2008b manages to reconcile some of the spectro- 



scopic ages of galactic clusters with CMD measurements when 
a hot star component is added to the fitting, to mimic the pres- 
ence of blue horizontal branch stars or blue stragglers. It is not 
straightforward, however, to predict how the HB morphology 
will impact the spectroscopic ages. The effective impact of the 
HB morphology on spectroscopic age is likely a non-trivial inter- 
play between the wavelength range fitted (in the sense that bluer 
regions will be more sensitive to hotter HB stars) and the exact 
morphology (or lack of) predicted by the underlying isochrone 
of the stellar population model, this morphology also being de- 
pendent on the age and metallicity of the mode lled population. 

Using the study by iGratton et al.l (1201 Oh on the galactic 
clusters HB morphologies, we searched for patterns of the HB 
morphologies that could correlate with the clusters with deviant 
ages, but could not find any. The only note of interest is that four 
of these clusters (NGC2808, NGC6388, NGC6441, NGC7078) 
have high values of (larger than 0.8). Nevertheless, cluster 
NGC 5927 (GC #9) also has a very high R' and its spectroscopic 
age matches the CMD range. 

Regarding the metallicities derived via spectral fitting, in 
Fig. [2] we compare our [Fe/H] values w ith those compiled by 
ICarretta et al] d2009l) (black triangles) and lSchiavon et al.l J2005I) 
(blue squares). The overall difference (this work - Schiavon et 
al.) is of -0.05 + 0.16 dex, while (this work - Carretta et al.) is - 
0. 15 + 0. 17 dex. Even though the agreement with Schiavon et al. 
scale is better on average, it was pointed out by the referee that 
there seems to be a slope between our results and Schiavon et 
al. scale (the relation crossing the line of equality). On the other 
hand, the offset between our results and Carretta et al. scale is 
nearly constant, thus the differences between our results and this 
scale might be explained by a zero-point offset. 

The scatter we obtain between our metallicities and the stel- 
lar analysis is compatible with the difference between Schiavon 
et al. and Carretta et al. scales (0. 10 + 0.20 dex), which are both 
based on high resolution stellar spectroscopy. It is a remarkable 
agreement between metallicities from integrated light at medium 
spectral resolutions and high-resolution stellar analysis. 

From the analysis of the galactic sample we conclude that re- 
garding the age determination, the method returned accurate res- 
ults except for six clusters: NGC 2808, NGC 5286, NGC 6121, 
NGC 6388, NGC 6441, NGC 7078. In the case of NGC 5286 
(and possibly NGC 6752, whose spectroscopic age marginally 
matches the CMD range), we suspect that contamination from 
foreground stars hampered the observations. NGC 2808 is a 
striking case of a cluster with triple main sequences, and a failure 
when comparing to SSP models is not disquieting. For the re- 
maining four clusters, or ~12% of our galactic sample, ages are 
underestimated by ~2-3 Gyr and we cannot provide a robust ex- 
planation for these differences. We confirmed that the presence 
of populations with CNONa anti-correlated abundances cannot 
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R' = N/y B /N' R0B , where N HB is the number of stars in the HB and 
N'/ios is the number of stars on the RGB brighter than V(HB)+1. 



-1.5 -1.0 -0.5 
[Fe/H] (this work) 

Figure 2. Top panel: Comparison of metallicities obtained with 
the spec tral fitting (this work) v ersus the values co mpiled in liter- 
ature b y TSchiavon et al.l (120051) (blue squares) and lCarretta et al.l 
d2009l) (black triangles) for the Galactic GCs. The error bar in 
the top left corner indicates the mean error values. Bottom panel: 
[Fe/H] residuals (ours minus literature values). 



explain the age differences. If this age difference is related to 
HB morphologies, as usually claimed in literature, its exact ef- 
fect remains to be better understood, given that not all GC with 
extended HB were affected. We reproduced the CMD ages of 
other GCs known to have blue components (such as NGC 5946 
and NGC 6284), without invoking addition al parameters in th e 
fitting (such as a free amount of hot stars in lKoleva et alJ 120081) . 
The mean difference (this work - literature) is of -0.8 Gyr with 
an r.m.s. of 1.7 Gyr, for the sub-sample with fitted ages inside 
the observational uncertainties; and -1.8 Gyr with a r.m.s. of 2.8 
Gyr for the whole sample. 

As for the metallicities, the method gives results in accord- 
ance with determinations from high resolution stellar spectro- 
scopy (R 30,000), with a r.m.s. of the same order of the r.m.s. 
between two different sets of high-resolution results as presen- 
ted in the previous paragraph. 

4.2. Stellar populations in M31 GCs 

We analysed the M3 1 sample with the same set up and procedure 
as in the galactic sample, and our derived ages and metallicities 
are presented in Table [3] We show in Fig. [3] the fitting of two 
spectra in our M31 sample, for illustration purposes. We obtain 
a large range of ages, from ~150 Myr (B322) to very old ages. 
In fact, three of the objects (B163, B393, and B39 8) were given 
"olde r than the Universe" ages (13.75 ± 0.11 Gvr: Ijarosik et al.l 
120111) . In these cases, x 1 maps of age distributions show a val- 
ley of low values starting around ~ 6 Gyr and almost flat with 
metallicity, indicating that the results are degenerate in age. The 
metallicities in our sample range from -2.2 dex to +0. 1 dex. 

Many of our GCs have been studied in literature, some by 
isochrones fitting to CMDs, or via spectral indices. A non- 
exhaustive list of results from literature is presented in Table [4] 
We compare our ages with those in the literature in Fig. [4] From 
this figure, we conclude there is, on average, general agreement 
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Figure 3. Spectral fitting of the GCs B163 (left-hand panel) and B380 (right-hand panel). The top panels show the spectra in black 
and the best fit in blue. The cyan lines are the multiplicative polynomials. The bottom panels are the residuals from the best fit, 
where the continuous green lines mark the 1-cr deviation as expected by the input S/N, and the dashed green lines mark the zero 
y-axis. This plot is best seen in colour (online version only). 



between our results and those in the literature, but with large 
dispersions. 

For the objects in common between this work and literature 
(see Table |4] and Fig 0s caption for details), we find a mean 
difference of 2.45 Gyr with a r.m.s. of 5.39 Gyr for photometry- 
based ages; and 0.06 Gyr with a r.m.s. of 4.70 Gyr for spectral 
indices ages. The agreement with metallicities (Pig. O is tighter: 
we determine a mean [Fe/H] difference of 0.21 dex (photometry- 
based, with a r.m.s. of 0.54 dex) and of 0. 10 dex with a r.m.s. of 
0.50 dex for those values determined through spectral indices. 
We also compared the age and metallicity mea surements for 
those M3 1 GCs in common with the studies of iBeaslev et al.l 
d2005l) and lPuzia et al.l d2005l) . Both groups have used the Lick- 
index method but differe nt SSP models. The me an age and 
[Fe/H] difference between [Beaslev et al.l d2005l) and iPuzia et al.l 
00051) results is of -0.60 Gyr (r.m.s. of 3.35 Gyr) and of -0.28 
dex (r.m.s. = 0.48 dex), respectively. 

We derive intermediate ag es (2-8 Gyr) for 7 GCs in our 
sample. ICaldwell et al.l d201 1 ) recently analysed a large sample 
of M31 clusters and they found that most of the intermedi- 
ate age clusters (as analysed by previous work in literature) 
were actually old metal-poor clu sters. The i r anal ysis is based 
on spectral indi ces and models by Schiavonl(l200"7l) . Using spec- 
tral fitting and IVazdekis et al.1 (l2010h SSP models, our results 
agree better with previous work which obtai n intermediate ages 
dBeaslev et al.ll200ll2005irPuzia et al.ll2005lf l 

We have three outer halo GC s in our sample th a t were pre- 
viously analysed in the study of lAlves-Brito et al 1 (I2009h . For 
MGC1, MGC5, and MGC 10 we find ages and metallicities of, 
respectively, (16.37, 9.99, and 11.60 Gyr) and (-1.36, -1.17, 
and -1.76 dex). Except for MGC1, for which we have found 
a higher age, there is a good agreement with the values of (7.10, 
10.00, and 12.60 Gyr) and (-1.3 7, -1.33, and -1.73 dex) repor- 
ted by lAlves-Brito et al.l d2009l) . respectively. MGC1 has also 
been recently investigated phot ometrically and spec troscopic - 
ally. Through the CMD analysis. iMackev et aT1(l2010l) conclude 
that MGC1 is old (12.5 to 12.7 Gyr) and metal-poor ([Fe/H] = 



4 But see as well lStrader et al] d2009h . who performed mass-to-light 
ratio analysis and favored old instead of intermediate ages for some of 
these GCs.. 



-2.3 dex). Spectrosco pically through ^-minimisation fittings, 
iFan & de Griisl d2012h also estimated an old age (13.3 Gyr) for 
MGC1 but average metallicities ranging from [Fe/H] = -2.06 to 
[Fe/H] = -1.76, depending on the SSP models used as well as 
on the stellar evolutio nary tracks employed. In addition, while 
IFan & de Griisl (1201 2|) ad opted a +0.42 dex enhan cement of a 
elements in their analysis, lAlves-Brito et al ] d2009l) have used a 
solar mixture, which could also account for the main differences 
in [Fe/H] between the different works. 

Photometric, spectral indices or spectral fitting methods dif- 
fer not only in their technical details, but also their underlying 
stellar population models differ in their choice of evolutionary 
tracks, the libraries of stellar s pectra and in their im plementation 
details, such as interpolations. ICoelho e t al. (2009), for example, 
made an analysis of the stellar population in M32 and has shown 
that using the same method with different SSP models result 
in different age and metallicity distributions; they conclude that 
choosing different SSP models from the literature might yield 
different results of age and metallicity at the quantitativ e level , 
even though a qualitative agreement is met. iDias et al ] (120101) 
also found evidence of this dependence of the results on the 
choice of SSP models when studying integrated spectra of GCs 
in the Small Magellanic Cloud. As we tested our method with 
GCs in the Milky Way, we believe that our techniques and mod- 
els are reliable, with the caveat that some objects (~12% from 
the galactic sample) might have their ages underestimated by 
~ 3 Gyr. 

IBeaslev et all (120051) ; IPuzia et al.l (120051) derive [c/Fe] val- 
ues as well as age and metallicities (in our Table |H we converte d 
their [Z/H] to [Fe/H] through equation 4 in lThomas et al]l2003h . 
Although there has been p rogress in measuring abundance pat- 
terns from spectral fitting dWalcher et al.ll2009D we have not at- 
tempted it in this work as the models currently available do not 
cover the parameter space needed. We intend to investigate this 
in a future work, with a larger sample. 

4.3. Comparing the two GC systems 

We show in Fig. [6] the age-[Fe/H] relation we obtain for our 
two samples, the Galactic GCs denoted by open squares and 
M31 GCs denoted by filled circles. It can be easily seen that 
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Table 3. Ages and metallicities for M31 GCs derived in this 
work. 



Cluster 


Age (Gyr) 


[Fe/H] 


B126 


2.41 +0.86 


-1.19 + 0.15 


B134 


13.72 ± 2.70 


-0.89 + 0.04 


B158 


12.44+ 1.24 


-0.74 + 0.02 


B163 


16.86 + 0.43 


-0.29 + 0.01 


B222 


1.16 + 0.12 


-0.28 + 0.03 


B225 


11.38 ± 1.74 


-0.44 + 0.03 


B234 


10.87 ± 1.10 


-0.73 ± 0.02 


B292 


4.09 + 0.47 


-1.54 + 0.04 


B301 


15.89 + 4.97 


-1.19 + 0.12 


B302 


11.23 ± 1.12 


-1.49 + 0.05 


B304 


8.53 ± 2.22 


-1.27 + 0.06 


B305 


0.98 + 0.10 


0.07 ± 0.02 


B307 


1.68 + 0.17 


0.02 + 0.02 


B310 


5.56 + 0.66 


-1.49 + 0.05 


B313 


12.67+ 1.27 


-0.83 + 0.02 


B314 


0.79 + 0.08 


-0.12 + 0.02 


B316 


1.46 + 0.26 


-0.00 + 0.13 


B321 


0.20 ± 0.02 


-0.07 + 0.06 


B322 


0.14 + 0.09 


-0.32 + 0.31 


B324 


1.00 + 0.10 


-0.13 + 0.09 


B327 


0.90 + 0.54 


-1.56 + 0.09 


B328 


2.58 + 0.81 


-1.60 + 0.05 


B331 


5.97 ± 1.34 


-0.64 + 0.07 


B337 


1.91+0.98 


-0.58 ± 0.32 


B347 


8.06+ 1.29 


-2.05 + 0.04 


B350 


8.70 ± 1.39 


-1.66 + 0.04 


B354 


11.43 ± 1.14 


-2.01 + 0.03 


B365 


9.01 ± 1.73 


-1.34 + 0.04 


B380 


0.58 + 0.06 


-0.06 + 0.03 


B383 


13.97 ± 1.37 


-0.57 + 0.03 


B393 


15.71 ± 1.58 


-1.00 + 0.02 


B398 


16.30+ 1.78 


-0.60 + 0.03 


B401 


8.49 + 0.85 


-2.22 + 0.04 


MGC1 


16.37 + 2.93 


-1.36 + 0.03 


MGC5 


9.99 ± 1.00 


-1.17 + 0.01 


MGC10 


11.60+ 1.16 


-1.76 + 0.01 


NB16 


2.90 + 0.56 


-1.09 + 0.11 


NB89 


12.14 + 2.42 


-0.77 ± 0.05 



the Galactic sample are consistent with old ages and a flat age- 
[Fe/H] relation, as expected from earlier studies on the GC 
system of our Galaxy. Using a sample of 54 Galactic GCs 
with high-quality CMDs and ages obtained photometrically, 
iFraix-Burnet et all (120091) found an age-metallicity relation for 
the halo GCs. According to their analysis the age-fFe/H] rela- 
tion goes from [Fe/H] = -2.0 dex 12 Gyr to [Fe/H] = -1.3 dex, 9 
Gyr ago. Nevertheless, their results are also consistent with a fast 
increase in metallicity (-1 .9 < [ Fe/H] < -1.4), at approxim ately 
10-11 Gyr ago. More recently. IForbes & Bridges! OOlOl) used 
n ew observations from the H ubble Space Telescope presented 
in lMarin-Franch et al.l d2009l) to compile a hig h-quality database 
of 93 Galactic GCs. With the larger sample, IForbes & Bridges! 
(1201 Oh showed that, in fact, the Galactic GCs can be divided 
into two groups. The first group is an old population of GCs, 
compatible with a rapid formation scenario of the Galactic halo, 
presenting a fiat age- [Fe/H] relation. The second group, how- 
ever, presents some younger objects displaying an age- [Fe/H] 
relation which is associated with disrupted dwarf galaxies. 



10.0 



a 
< 




1.0 10.0 
Age. Gyr (this work) 

Figure 4. Comparison between our derived ages and those pub- 
lished in literature for GCs in M31. Each coloured-symbol 
correspond to a different literature work, fol lowing Table 
[4] (a) blue diam ond correspond to resul ts by iBeaslev et al.l 
120051) ( models bvlBruzual & Charlolll2003l): (b) red upward tri- 
angle s to IBeaslev et al.l d2005l) (models by iThomas et alJl2003l 
12004 : (c) green squares to lPuziaetal.1 (120051) : (d) brown 
crosses to Fan et al] (20061) : (e) magenta downwa rd triangles 



to Beaslev et al. (I2004h~ (f) cyan o pen circles to | Wang et al.l 
d2010l) , and: (g) grey filled circles to lAlves-Brito et alj d2009lh 
As labeled in Table [4] references (a-c) and (e) present the res- 
ults based on photometric measurements, while references (d), 
(f) and (g) above present the results based on spectral indices. 
This plot is best seen in colour (online version). 



However, while the Galactic GCs are predominantly old, 
M31 shows not only young objects (< 1 Gyr), but also a 
number of interm ediate age and old GCs (Be aslev et al.ll200"5l : 
iPuzia etal 1120051) . These findings are also confirmed in our ana- 
lysis as shown in Table [3] From Fig. [6] (filled symbols), we see 
that, within the uncertainties, both old GCs in the Galaxy and in 
M31 present similar flat age-[Fe/H] relations. 

By contrast, we see that M31's intermediate/young popula- 
tion (age < 8 Gyr) follows a roughly linear age-[Fe/H] relation 
in the sense that the young GCs are also more metal-rich. These 
remarkable differences in age- [Fe/H] between the Galactic and 
M3 1 GC systems suggest that while the latter has likely experi- 
enced a recent, active merger history, the former has instead ex- 
per ienced a quiet star forma tion history over the last 10 Gyr (but 
see lForbes & Bridgesl2010l for a recent discussion on the topic). 
We notice, however, that a larger sample of young GCs in M3 1 
has to be targeted to confirm (or not) the age- [Fe/H] relation we 
recover in this work. 

The relation between the GC system's stellar population and 
kinematic properties can s hed light into forma tion of structures 
in a galaxy. For example, iPerrett et al.l d2002l) investigated the 
kinematic properties of several hundred GCs in M3 1 . They con- 
cluded that the metal-rich GCs present a centrally concentrated 
spatial distribution with a high rotation amplitude, consistent 
with a bulge population. Whereas the metal-poor GCs tend to be 
less spatially concentrated and were also fo und to have a strong 
rotation signature (see as well discussions in lHuchra et alJlT991c 
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-2.5 -2.0 -1.5 -1.0 -0.5 0.0 
[Fe/H] (this work) 
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Figure 5. Comparison between our derived metallicities [Fe/H] 
and those published in literature for GCs in M3 1 . The colour and 
symbol codes are the same as in Fig. [4] This plot is best seen in 
colour (online version). 
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Figure 6. Age-metallicity relation obtained for M31 (filled 
circles) and Galactic GCs (open squares). Refer to the text of 
the paper for more details. 



lBekkil2010tlLee et al.l2008tlMorrison et alj201 ll and references 
therein). 

Also, the dis(similarities) between the metallicity distribu- 
tion function (MDF) of the Galactic and M3 1 GC systems have 
been widely debated in literature. In Fig. [7] we present the MDF 
we recover from our own analysis for both systems. To our 
knowledge, this is the first direct comparison of [Fe/H] distri- 
butions in both Galactic and M31 GC systems employing the 
same methods and techniques, regardless of the GCs position 
in the galaxy. In the top panel, we see that the MDF obtained 
from our integrated metallicities for the Galactic GCs is cl early 
bimodal, a result that goes back to the seminal paper by Zinn 



(1 19851) . who was the first to propose the exist ence of two sub- 
populations in the Galactic GC system (see also lBica et al .120061 
for a recent discussio n). The KMM mixture modelling algorithm 
(I Ashman et al.ll 19941) suggests statistically convincing evidence 
of bimodality in the Galactic GC system (at better than the 99% 
confidence level). 

For M31, however, the shape and distribution of the MDF is 
still controversial. While some authors propose that the MDF of 
M3 1 GCs presents two sub-populations - one with a metallicity 
peak at [Fe/H] = -1.57 that is associated to the galaxy's halo, 
and other one peaking at [Fe/H] = -0.6 1, which is structurally 
associated to the ga l axy's bulge (see , e.g.. [Ashman & Birdll993i: 
iBarmbv et al.ll2000t iFan et al j2008T) - other authors suggest that 
the bimodality is not present at all or it is weakly detected (e.g. 
ICaldwell et alllloTlh . 

From our own analysis for both Galactic and M3 1 GCs using 
the same methods and techniques, we obtain that both systems, 
regardless the age differences, cover approximately the same 
[Fe/H] range. Within our limited sample, the old M31 GCs ex- 
tend from [Fe/H] = -2.2 up to [Fe/H] = -0.3, while the younger 
population goes from [Fe/H] = -2.2 up to 0. 1 dexQ Furthermore, 
the KMM algorithm applied to the old GCs in M3 1 does not sup- 
port bimodality (the probability is less than 77%), with the num- 
ber of metal-poor and metal-rich objects being almost the same 
in the galaxy. We note, however, that M3 1 's GC system is more 
than a factor of 2 larger than the Milky Way's. In addition, our 
sample is biased not only by the number of objects studied but 
also by their position in the galaxy. While our sample is biased 
towards disk/bulge objects, more GCs in the halo of M3 1 need to 
be targeted to better understand M3 1 's MDF shape and properly 
compare it with that of the Galaxy. 



5. Summary and Conclusions 

Spectroscopic ages and metallicities were derive d for a sample 
of 38 G Cs in M3 1 , drawn from the o bservations of lBeaslev et alJ 
(12001) andlXl ves-Brito et al.l d2009l) . These parameters were ob- 
tai ned by fitting the obse rved integrated spectra to SSP models 
by Vazdekis et alJ jl2010h using the spectral fitting code ULySS 
dKoleva et al.ll2009l) . To our knowledge, this is the first time that 
full spectrum fitting is used in deriving stellar population para- 
meters in M31 GCs. 

We tested the reliability of our analysis by fitting the in- 
tegrate d spectra of Galactic GCs presented in ISchiavon et al.l 
(120051) . In six cases, out of 34 objects for which we obtained 
CMD ages from the literature, the spectroscopic ages do not 
match the ages drawn from CMD analysis. In the case of 
NGC 5286 (and possibly NGC 6752), we suspect that contam- 
ination from foreground stars hampered the observations. This 
is unlikely to be an issue for extra-galactic clusters. NGC 2808 
is a striking case of a cluster with triple main sequences and 
complex HB morphology, and would deserve a more detailed 
modelling than SSP fitting. For the remaining four clusters, ages 
are underestimated by ~2-3 Gyr and we cannot provide a robust 
explanation for these differences. We did not find evidence of a 
correlation with either contamination from CNONa abundance 
variations or a specific HB morphology. 

The spectroscopic integrated metallicities der ived with spec- 
tral fitting were compared to the compilations bv lSchiavon et al.l 



Note that the lowest metallicity in the SSP models grid is [Fe/H] = 
-2.3. Therefore, if there is tail in the MDF extending to lower metalli- 
cities, we would not be able to see it in this analysis. 
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Figure 7. The metallicity distribution function for (a) Galactic 
GCs, (b) old (ages > 8 Gyr) M31 GCs, and (c) all (old+young) 
M31 GCs as analysed in this work in bins of 0.40 dex wide. 



(120051) and lCarretta et al.1 (120091) . Our results agree well with the 
higher resolution stellar analysis. 

As for M31, we obtain a large range of ages (from ~150 Myr 
to the age of the Universe) and metallicities (-2.2 < [Fe/H] < 
+0.1). We confirm previous results in the literature that find 
young globular clusters in M31, in contrast with the globular 
cluster system in our own Milky Way. We find an age-metallicity 
relation that can be described as having two components: an 
old population with a flat age- [Fe/H] relation, possibly associ- 
ated with the halo and/or bulge in analogy to what is seen in 
the Milky Way, and a second one with a roughly linear relation 
between age and metallicity (higher metallicities corresponding 
to younger ages). 

Finally, our analyses do not support [Fe/H] bimodality in the 
M3 1 globular cluster system, but further investigation with a lar- 
ger and unbiased sample is necessary to better interpret the star 
formation history of M3 1 when compared to our own Galaxy. 

Acknowledgements. EC and PC are thankful to Mina Koleva for her con- 
stant help on using the ULySS package. PC acknowledges financial support 
from FAPESP project 2008/58406-4. AB acknowledges support from Australian 
Research Council (Super Science Fellowship, FS1 10200016). JB acknowledges 
support from NSF grants AST-121 1995 and AST-1 109878. 



References 

Alves-Brito, A., Forbes, D. A., Mendel, J. T., Hau, G. K. T., & Murphy, M. T. 

2009, MNRAS, 395, L34 
Alves-Brito, A., Hau, G. K. T., Forbes, D. A., et al. 2011, MNRAS, 417, 1823 
Ashman, K. M. & Bird, C. M. 1993, AJ, 106, 2281 



Ashman, K. M., Bird, C. M., & Zepf, S. E. 1994, AJ, 108, 2348 
Barmby, P., Huchra, J. P., Brodie, J. P., et al. 2000, AJ, 1 19, 727 
Beasley, M. A., Brodie, J. P., Strader, J., et al. 2004, AJ, 128, 1623 
Beasley, M. A., Brodie, J. P., Strader, J., et al. 2005, AJ, 129, 1412 
Beaulieu, S. F, Gilmore, G, Elson, R. A. W., et al. 2001, AJ, 121, 2618 
Bekki, K. 2010, MNRAS, 401, L58 

Bica, E., Bonatto, C, Barbuy, B., & Ortolani, S. 2006, A&A, 450, 105 
Brodie, J. P. & Huchra, J. P. 1990, ApJ, 362, 503 
Bruzual, G. & Chariot, S. 2003, MNRAS, 344, 1000 

Buonanno, R., Corsi, C. E., Pulone, L., Fusi Pecci, F, & Bellazzini, M. 1998, 
A&A, 333, 505 

Caldwell, N., Schiavon, R., Morrison, H, Rose, J. A., & Harding, P. 201 1, AJ, 
141, 61 

Carretta, E., Bragaglia, A., Gratton, R., D'Orazi, V., & Lucatello, S. 2009, A&A, 
508, 695 

Carretta, E. & Gratton, R. G. 1997, A&AS, 121, 95 

Cenarro, A. J., Peletier, R. F, Sanchez-Blazquez, P., et al. 2007, MNRAS, 374, 
664 

Chaboyer, B., Demarque, P., Kernan, P. J., & Krauss, L. M. 1998, ApJ, 494, 96 
Cid Fernandes, R. & Gonzalez Delgado, R. M. 2010, MNRAS, 403, 780 
Coelho, P., Bruzual, G., Chariot, S., et al. 2007, MNRAS, 382, 498 
Coelho, P., Mendes de Oliveira, C, & Cid Fernandes, R. 2009, MNRAS, 396, 
624 

Coelho, P., Percival, S., & Salads, M. 2012, ArXiv e-prints 

Coelho, P., Percival, S. M., & Salads, M. 2011, ApJ, 734, 72 

Conroy, C. & Gunn, J. E. 2010, ApJ, 712, 833 

De Angeli, F, Piotto, G., Cassisi, S., et al. 2005, AJ, 130, 1 16 

de Freitas Pacheco, J. A. & Barbuy, B. 1995, A&A, 302, 718 

Delgado, R. M. G., Cervino, M., Martins, L. P., Leitherer, C, & Hauschildt, P. H. 

2005, MNRAS, 357, 945 
Dias, B., Coelho, P., Barbuy, B., Kerber, L., & Idiart, T. 2010, A&A, 520, A85+ 
Fan, Z. & de Grijs, R. 2012, MNRAS, 424, 2009 

Fan, Z., Ma, J., de Grijs, R., Yang, Y., & Zhou, X. 2006, MNRAS, 371, 1648 
Fan, Z., Ma, J., de Grijs, R., & Zhou, X. 2008, MNRAS, 385, 1973 
Forbes, D. A. & Bridges, T. 2010, MNRAS, 404, 1203 
Fraix-Burnet, D., Davoust, E., & Charbonnel, C. 2009, MNRAS, 398, 1706 
Gallart, C, Zoccali, M., & Aparicio, A. 2005, ARA&A, 43, 387 
Galleti, S., Federici, L., Bellazzini, M., Fusi Pecci, F, & Macrina, S. 2004, A&A, 
416, 917 

Girardi, L., Bressan, A., Bertelli, G., & Chiosi, C. 2000, A&AS, 141, 371 

Gratton, R. G., Carretta, E., & Bragaglia, A. 2012, A&A Rev., 20, 50 

Gratton, R. G., Carretta, E., Bragaglia, A., Lucatello, S., & D'Orazi, V. 2010, 

A&A, 517, A81 
Holland, S. 1998, AJ, 115, 1916 

Huchra, J. P., Brodie, J. P., & Kent, S. M. 1991, ApJ, 370, 495 

Jarosik, N, Bennett, C. L., Dunkley, J., et al. 2011, ApJS, 192, 14 

Kent, S. M., Huchra, J. P., & Stauffer, J. 1989, AJ, 98, 2080 

Koleva, M., Prugniel, P., Bouchard, A., & Wu, Y. 2009, A&A, 501, 1269 

Koleva, M., Prugniel, P., Ocvirk, P., Le Borgne, D., & Soubiran, C. 2008, 

MNRAS, 385, 1998 
Kraft, R. P. & Ivans, 1. 1. 2003, PASP, 1 15, 143 

Le Borgne, D., Rocca-Volmerange, B., Prugniel, P., et al. 2004, A&A, 425, 881 

Lee, H.-c, Worthey, G., Dotter, A., et al. 2009, ApJ, 694, 902 

Lee, H.-c, Yoon, S.-J., & Lee, Y.-W. 2000, AJ, 120, 998 

Lee, M. G., Hwang, H. S., Kim, S. C, et al. 2008, ApJ, 674, 886 

Ma, J., Yang, Y, Burstein, D., et al. 2007, ApJ, 659, 359 

Mackey, A. D., Huxor, A. P., Ferguson, A. M. N., et al. 2010, ApJ, 717, LI 1 

Maraston, C. 2005, MNRAS, 362, 799 

Marin-Franch, A., Aparicio, A., Piotto, G, et al. 2009, ApJ, 694, 1498 
Meissner, F. & Weiss, A. 2006, A&A, 456, 1085 

Mendel, J. T., Proctor, R. N, & Forbes, D. A. 2007, MNRAS, 379, 1618 
Milone, A. D. C, Sansom, A. E., & Sanchez-Blazquez, P. 2011, MNRAS, 414, 
1227 

Momany, Y, Ortolani, S., Held, E. V., et al. 2003, A&A, 402, 607 
Morrison, H, Caldwell, N, Schiavon, R. P., et al. 2011, ApJ, 726, L9 
Ocvirk, P. 2010, ApJ, 709, 88 

Percival, S. M., Salaris, M., Cassisi, S., & Pietrinferni, A. 2009, ApJ, 690, 427 
Perrett, K. M., Bridges, T. J., Hanes, D. A., et al. 2002, AJ, 123, 2490 
Piotto, G., Bedin, L. R., Anderson, J., et al. 2007, ApJ, 661, L53 
Proctor, R. N, Forbes, D. A., & Beasley, M. A. 2004, MNRAS, 355, 1327 
Puzia, T. H, Perrett, K. M., & Bridges, T. J. 2005, A&A, 434, 909 
Rutledge, G. A., Hesser, J. E., Stetson, P. B., et al. 1997, PASP, 109, 883 
Sanchez-Blazquez, P., Peletier, R. F, Jimenez- Vicente, J., et al. 2006, MNRAS, 
371, 703 

Schiavon, R. P. 2007, ApJS, 171, 146 

Schiavon, R. P., Rose, J. A., Courteau, S., & MacArthur, L. A. 2004, ApJ, 608, 
L33 



9 



Cezario et al.: Spectral fitting of globular clusters 



Schiavon, R. P., Rose, J. A., Courteau, S., & MacArthur, L. A. 2005, ApJS, 169, 
163 

Strader, J., Smith, G. H., Larsen, S., Brodie, J. P., & Huchra, J. P. 2009, AJ, 138, 
547 

Thomas, D., Maraston, C, & Bender, R. 2003, MNRAS, 339, 897 
Thomas, D., Maraston, C, & Korn, A. 2004, MNRAS, 351, L19 
Usher, C, Forbes, D. A., Brodie, J. P., et al. 2012, MNRAS, 426, 1475 
Vazdekis, A., Sanchez-Blazquez, P., Falcon-Barroso, J., et al. 2010, MNRAS, 
404, 1639 

Walcher, C. J., Coelho, P., Gallazzi, A., & Chariot, S. 2009, MNRAS, 398, L44 

Wang, S., Fan, Z., Ma, J., de Grijs, R., & Zhou, X. 2010, AJ, 139, 1438 

Wolf, M. J., Drory, N, Gebhardt, K., & Hill, G. J. 2007, ApJ, 655, 179 

Worthey, G., Faber, S. M., Gonzalez, J. J., & Burstein, D. 1994, ApJS, 94, 687 

Zinn, R. 1985, ApJ, 293, 424 

Zinn, R. & West, M. J. 1984, ApJS, 55, 45 

Zoccali, M., Renzini, A., Ortolani, S., et al. 2001, ApJ, 553, 733 



10 



Cezario et al.: Spectral fitting of globular clusters 
Table 2. Ages and metallicities for Galactic GCs. 
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Table 4. Ages and metallicities from literature, for our sample of M31 GCs. 
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OpcULld.1 lllUICCa 


a 




A 90 
D.ZU 


-1-110 
it 1 . 1 u 


1 97 _>_ SS 
- 1 .Z / it U.J J 


Spectral indices 


U 




A 1 
L t. 1 U 


-1- 1 70 
± J. /U 


SO -1- 9A 
-U.JU it U.ZO 


Spectral indices 


c 








n 7A -1- n js 
-U. /o ± U.Zj 


Photometry 


f 
1 


R^nA 


1 9 QO 
IZ.VU 


-1- A 70 


1 ss -1- n as 

-1 . JJ ± U.OJ 


Spectral indices 


a 




1 1 QO 

1 i .yu 


-1- ^ QO 
± J.VU 


1 /is -1- n 

-1 .40 ± U.4J 


Spectral indices 


D 




1 ^ 90 


It J. 1U 


1 1 A -1- 4A 


Vv / * 1 1" ' 1 1 i n/11 r*pc 
opCULIal lllUlCCa 


c 








1 Li -1- n jo 
-l.JZ ± u.zz 


Photometry 


1 


R^ns 


AO 
U.4U 


-1- n 1 n 

± U. 1U 


n on -1- n ai 

-u.vu ± U.Ol 


Photometry 


f 

1 


rs JU / 


1 A1 
1.D1 


-1- n 1 n 

± U. 1U 


n a 1 -1- n ^a 

-U.41 ± U.JO 


Photometry 


f 
I 


R^ 1 


1 Q SO 


+ 910 


9 OA -1- s 1 

-Z,.UU It U.J 1 


v» / * 1 1" ' 1 1 i n/ii r*pc 
OpCULIal lllUlCCa 


a 




1 1 /in 

1 1 .4U 


-u "3 SO 

± j.jU 


1 AO -1-0/17 
-1.0V ± U.4/ 


Spectral indices 


D 




1 n sn 

1U.0U 


j. q in 

it J. 1U 


1 9A -1- n S9 

- 1 .ZH- it U.J Z 


Spectral indices 


c 








1 Al -a- n 98 

- 1 .H J it U.ZO 


Photometry 


f 

1 


R31 3 


s sn 

o.oU 


■ 1 on 
± 1 .zu 


1 ns -1- n ^s 

- 1 .UJ ± U.J J 


Spectral indices 


a 




1 1 7n 


-1- n on 
± u.yu 


n so -1- n sn 
-u.ov ± u.ju 


Spectral indices 


D 




1 1 9n 
1 1 .zu 


-1- 1 9n 
± 1 .zu 


n 70 -1- n ^9 

-U. Iy± U. JZ 


Spectral indices 


C 




7 98 
/ .ZO 


± u. /u 


1 no -1- n 1 n 

- 1 .U7 ± U. 1U 


Photometry 


f 
1 


R3 1 A 


n sn 


-1- n ^.n 
± u.ou 


n ^s -1- n 9s 
U.JJ ± U.Zj 


Spectral indices 


c 




1 nn 

1 .UU 


-1- n 1 n 
± u. 1U 


n ^s -1- n 99 
-U.JJ ± u.zz 


Spectral indices 


c 


R3 1 A 
ID J 10 


1 nA 
i.uo 


-1- n 1 n 

± U. 1U 


1 A7 -1- n 9^ 

-1 .4 / ± U.ZJ 


Photometry 


f 
I 


R^9 1 


n ^n 

U.JU 


+ n ^n 

± U.JU 


n nn -1- n An 

U.UU ± u.ou 


Spectral indices 


c 




1 nn 
1 .uu 


-1- n 1 n 

± U. 1U 


1 os _i_ n ^n 
-1,70 ± U.JU 


Spectral indices 


c 


R399 
D JZZ 


n 1 n 

U. 1U 


-1- n sn 

± U.JU 


n 9n -1- n 9n 
-u.zu ± u.zu 


Spectral indices 


c 




Z. jU 


-u n 7n 
± u. /u 


1 QQ -1- n A L 
-L.yy ± U.4J 


Spectral indices 


c 


R39A 
D JZ4 


n on 
u.vu 


-1- n 9n 
± u.zu 


n ns -1- n An 

-U.UJ ± U.4U 


Spectral indices 


c 




1 nn 

l.UU 


-1- n 1 n 

± U. 1U 


n ss -1- n 9n 
-U.Jo ± u.zu 


Spectral indices 


c 


R^97 


n 1 n 

U. 1U 


-1- n on 

± U.VU 


n ^n -1- n 7S 

U.JU ± u. / J 


Spectral indices 


c 




5.40 


+ 1.40 


-1.97 + 0.34 


^lnpptr!nl inHipp*; 




B328 


16.80 


±5.20 


-2.13 ±0.35 


Spectral indices 


a 




13.30 


± 1.30 


-1.69 ±0.35 


Spectral indices 


b 




11.00 


±2.40 


-2.29 ± 0.42 


Spectral indices 


c 


B337 


2.60 


± 1.90 


-1.18 ± 0.19 


Spectral indices 


a 




6.60 


±3.20 


-1.19 ±0.25 


Spectral indices 


b 




4.90 


±2.90 


-0.59 ±0.11 


Spectral indices 


c 




2.03 


±0.10 


-1.09 ±0.32 


Photometry 


f 



Continued on next page 
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Table 4. continued. 



Cezario et al.: Spectral fitting of globular clusters 



Cluster 




Age 
(Gyr) 




[Fe/H] 

(dex) 


Method 


Reference 


B347 


8.70 


± 5.70 


-2.49 


± 0.21 


Spectral indices 


a 




9.90 


±4.00 


-1.93 


± 0.51 


Spectral indices 


b 




10.20 


± 2.40 


-2.61 


± 0.37 


Spectral indices 


c 




2.53 


± 0.15 


-1.71 


± 0.03 


Photometry 


f 


B350 


9.80 


± 2.50 


-1.98 


± 0.49 


Spectral indices 


a 




12.40 


±4.70 


-1.71 


± 0.33 


Spectral indices 


b 




9.30 


± 2.30 


-1.81 


± 0.37 


Spectral indices 


c 




1.99 


± 0.10 


-1.47 


± 0.17 


Photometry 


f 


B354 


5.24 


± 0.65 


-1.46 


± 0.38 


Photometry 


f 


B365 


9.20 


± 3.10 


-1.59 


± 0.44 


Spectral indices 


a 




10.40 


± 2.20 


-1.42 


± 0.51 


Spectral indices 


b 




6.60 


± 3.10 


-1.03 


± 0.20 


Spectral indices 


c 




1.73 


± 0.10 


-1.78 


± 0.19 


Photometry 


f 


B380 


0.45 


± 0.13 


0.15 


± 0.10 


Spectral indices 


c 




1.00 


± 0.10 


-1.69 


± 0.13 


Spectral indices 


c 


B383 


10.30 


± 2.20 


-0.68 


± 0.18 


Spectral indices 


a 




10.60 


± 1.30 


-0.68 


± 0.31 


Spectral indices 


b 




11.30 


± 2.40 


-0.59 


± 0.08 


Spectral indices 


c 




13.99 


± 1.05 


-0.48 


± 0.20 


Photometry 


f 


B393 


10.00 


± 1.80 


0.00 


± 0.00 


Spectral indices 






11.90 


± 1.30 


-0.80 


± 0.47 


Spectral indices 


b 




11.20 


± 1.40 


-0.60 


± 0.14 


Spectral indices 


c 




6.76 


± 1.10 


-1.41 


± 0.05 


Photometry 


f 


B398 


12.10 


± 1.50 


-0.71 


± 0.41 


Spectral indices 






14.70 


± 3.70 


-0.52 


± 0.28 


Spectral indices 


b 




15.00 


± 2.40 


-0.54 


± 0.10 


Spectral indices 


c 


B401 


9.20 


± 5.80 


-2.49 


± 0.49 


Spectral indices 


a 




15.00 


± 5.10 


-2.38 


± 0.24 


Spectral indices 


b 




10.20 


± 2.40 


-2.38 


± 0.32 


Spectral indices 


c 




3.49 


± 0.40 


-1.75 


± 0.29 


Photometry 


f 


MGC1 


7.10 


± 3.00 


-1.37 


± 0.15 


Spectral indices 




MCGC5 


10.00 


± 3.00 


-1.33 


±0.12 


Spectral indices 


g 


MCGC10 


12.60 


±3.00 


-1.73 


±0.20 


Spectral indices 


g 


NB16 


2.00 


± 1.40 


-1.40 


±0.18 


Spectral indices 


a 




4.60 


±0.80 


-1.33 


±0.24 


Spectral indices 


b 




7.20 


±3.50 


-1.11 


±0.15 


Spectral indices 


c 


NB89 


8.30 


±2.60 


-0.78 


±0.14 


Spectral indices 


a 




14.40 


± 2.70 


-0.84 


±0.15 


Spectral indices 


b 




6.80 


± 1.70 


-0.58 


±0.08 


Spectral indices 


c 



( models by 



(d) iFan et al 



Bruzual & Charlotll2003h; ( b) IBeaslev et al] d2005h (mod els by 
j2006T ): (e) IBeaslev et al] d2004l) ; (f) IWang et al] d2010h . and; 



References for column 5 : (a) IBeaslev et al] d2005 
iThomas et ai]|2003L 120041); (c) iPuzia et al] d2005l) ; 

fg ) lAlves-Brito et al] 12009). Regarding the references (a), (b) and (c), we converted their [Z/H] to [Fe/H] through equation 4 
in lTholr!aTeTaU "i r2003l) . Notes: We did not find in the listed literature the parameters for two clusters in our sample, B302 and 
B331. 



